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Abstract 


The main goal of this work is to develop and analyze an accurate trun- 
cated stochastic Runge-Kutta (TSRK2) method to obtain strong numeri- 
cal solutions of nonlinear one-dimensional stochastic differential equations 
(SDEs) with continuous Hélder diffusion coefficients. We will establish 
the strong L!-convergence theory to the TSRK2 method under the local 
Lipschitz condition plus the one-sided Lipschitz condition for the drift co- 
efficient and the continuous Hélder condition for the diffusion coefficient 
at a time T and over a finite time interval [0,T], respectively. We show 
that the new method can achieve the optimal convergence order at a finite 
time T compared to the classical Euler-Maruyama method. Finally, nu- 
merical examples are given to support the theoretical results and illustrate 
the validity of the method. 
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1 Introduction 


Consider the scalar stochastic differential equation of It6 type 


au())dg + W(a(E))dW(E), — € € [0,7], (1) 
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where T > 0 and a,b: R > R are Borel measurable functions and W(€) 
is a one-dimensional Wiener process. We refer the reader to [5, 16, 24, 26] 
for an overview of stochastic differential equations (SDEs) and their applica- 
tions. We assume that the drift coefficient a satisfies the one-sided Lipschitz 
and local Lipschitz conditions and that the diffusion coefficient b satisfies the 
Holder continuity condition. Such applications exist in finance, for example, 
the Ait-Sahalia-type interest rate model [1, 4] and the Cox—Ingersoll—Ross 
model [6], and in biology, for example, the stochastic SIS epidemic model 
[7] and laser emission models in chemical physics [16]. Since these equations 
can be very complicated and analytical solutions are not always possible, 
numerical methods have become an efficient tool for computing approximate 
solutions for SDEs. Many numerical methods have been developed for SDEs 
under the global Lipschitz and linear growth conditions; see, for example, 
(24, 16, 29]. In particular, derivative-free stochastic Runge-Kutta (SRK) 
methods for strong approximations have been proposed [3, 27, 8]. How- 
ever, Hutzenthaler, Jentzen, and Kloeden [13] have shown that the classical 
Euler-Maruyama (EM) and Milstein methods do not converge strongly to 
the solution of (1) when the global Lipschitz and linear growth conditions of 
the drift or diffusion coefficients are perturbed. 


In recent years, an increasing number of numerical methods have been 
developed for solving nonlinear SDEs without global Lipschitz conditions. 
These methods include implicit methods [25, 23, 9], the tamed numerical 
methods, the first of which was presented in Hutzenthaler et al. [14, 28, 12], 
and the stopped EM method [19]. However, the use of implicit or drift- 
implicit numerical methods requires the solution of a nonlinear algebraic 
equation at each time step and thus can be very inefficient. Moreover, the 
tamed methods may lead to inaccurate results due to the perturbation of the 
flow by changing the drift and diffusion coefficients even at moderately small 
step sizes [30]. Recently, Mao [21] proposed the truncated EM method for a 
strong approximation of the nonlinear SDEs under the local Lipschitz con- 
dition and the Khasminskii-type condition. After that, the [’-convergence 
rates and stability properties of the truncated EM method have been stud- 
ied by some researchers [22, 11]. Yang et al. [32] investigated the strong 
convergence of the truncated EM method for one-dimensional SDEs with 
superlinearly growing drifts and the Hélder continuous diffusion conditions. 
Then, several new techniques of the partially truncated EM method were 
proposed in [33] to determine the optimal convergence rate. The authors 
also investigated the stability of these methods. 


Despite the strength of explicit methods in terms of computational cost, 
there is still a drawback. As mentioned in [16, 2], explicit methods may have 
to use very small step sizes if the SDEs to be solved are stiff. Although under 
the classical global Lipschitz condition, there are some classes of explicit 
methods with extended stability regions that are well suited for solving stiff 
problems, especially those whose eigenvalues are close to the negative real 
axis. For example, Komori and Burrage [17] have developed strong first- 
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order SROCK methods for It6 and Stratonovich SDEs. However, they are 
still far from accurate numerical algorithms suitable for highly nonlinear stiff 
SDEs with explicit methods. 


In this paper, we will bring all these ideas together. Based on Mao’s 
truncation strategies [21, 22, 18], we derive a stochastic two-stage truncated 
Runge-Kutta method for nonlinear SDEs with the superlinearly growing drift 
coefficient and the continuous Hélder diffusion coefficient. The proposed 
scheme is explicit and includes some free parameters that can extend the 
accuracy of the results and stability regions. To the best of the author’s 
knowledge, truncated SRK is the only truncated Runge-Kutta method with 
the relevant results for a strong approximation of solutions of SDEs with 
Holder diffusion coefficients. We will study the strong convergence of the 
proposed method under the local Lipschitz condition plus the one-sided Lip- 
schitz condition for the drift term and the continuous Hélder condition for the 
diffusion term at a time T and over a finite time interval [0, T], respectively. 
We show that the new method can achieve the optimal order of convergence 
compared to the classical EM method at a finite time T without any restric- 
tion on the step size. To show the effectiveness of our methods, we simulate 
some stiff SDEs with several Hélder parameters a € [0, 3). 


The rest of the paper is organized as follows. In Section 2, we describe 
some relevant assumptions that must be satisfied for the drift and diffusion 
coefficients, as well as results for solving the original SDEs. In Section 3, we 
first present mathematical notations and preliminary results for truncated 
methods. We then develop the two-stage truncated stochastic Runge-Kutta 
(TSRK2) method, which is the main goal of the paper. In this section, we 
present convergence results of the new method at time T and give several 
technical lemmas. In Section 4, we study the convergence rate over a finite 
interval [0,7]. Numerical results and conclusions are given in Sections 5 and 
6, respectively. 


2 Mathematical preliminaries 


Let (Q,F7,P) be a complete probability space with right continuous and 
increasing filtration {Fe}o<e<r, where Fo contains all P-null sets. Let 
{W(€)}o<e<r be a one-dimensional standard {F¢}o<e<r-adapted Wiener 
process on the probability space. For 21, z2 € R, we use 21 V 22 = max{2z, za} 
and 21 A z2 = min{z1, z2}. If Bis a set, then its indicator function is denoted 
by Ig, namely, Ip(z) = 1 if z € B and 0 otherwise. For any fixed p € [1, oo), 
we frequently make use of the Young inequality 


5 - _ 
oes maa + Eee zpla for all 21,2 € R™, (2) 
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for any 6 > 0 and q € [p,oo). To construct the new method, we now make 
assumptions about a and b. 


Assumption 1. Suppose that there exist real positive constants Ky, and p 
such that 


|a(z1) — a(z2)|? < Ki(1 + Jarl? + |22l")|z1 — 22l?, (3) 
for all 21, z2 € R%. 


From (3) we can conclude that the drift coefficient a satisfies the local 
Lipschitz condition: For any u > 0, there exists K, > 0 such that 


|a(z1) — a(z2)| < Kuler — 22], (4) 
for all 21, 22 € R with |z| V |z| <u. 


Assumption 2. We assume that the drift coefficient a satisfies the one- 
sided Lipschitz condition and the diffusion coefficient b satisfies the Hélder 
continuity condition: Real constants H,,H2 € Rt and 0 < a < 0.5 exist 
such that 


(z1 — 22)(a(z1) — a(z2)) < Mila — 29’, (5) 
|b(z1) — b(z2)| < Halex — 20/2 *%, (6) 


for all z1,z2 € R. 
Clearly, from Assumption 2, it is easy to verify that 


za(z) < Hy|lz2|? + |z||a(0)| < Mi(1 + |2|*), (7) 
|b(z)| < He|z|°F* + |b(0)| < Mo(1 + |2\), (8) 


for all z € R, where M, = 0.5(|a(0)|? V (2H; + 1]) and Mz = H2(0.5+ a) V 
[H2(0.5 — a) + |b(0)|]. The relation (8) shows that the coefficient function b 
satisfies the linear growth condition. 

In this work, we use C for the generic positive real constant that depends 
on p, T, a, Xo, and so on but are independent of the time step size A and R, 
and whose values can change between occurrences. 


Remark 1. Under Assumption 2, for all p € (2,00), there is a positive 
constant C' such that 


—1 
za(z) + lo P <C(1+|z|?), for all zeER. (9) 
To construct the numerical method, we first estimate the growth rate of a 


as follows. We choose a strictly increasing continuous function v : Rt > Rt 
such that v(r) + co as r + oo and 
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sup laa) = a(z2)/ V sup |a(z)| < v(r), forallr >1. (10) 


O<|zi|Vizal<r [21 — 2al lz|<r 


We see that v~! : [v(0),00) > (0,00) as the inverse function of v is also a 
strictly increasing continuous function. From Assumption 1, we observe that 


sup jazi) — a(za)] 4 J/Ki(1 + V2re/?), 


0<|z1|V|z2|<r |21 = z2| 
sup |a(z)| < Ja(0)| + V/Kir(1 + r?/?), (11) 
lzl<r 


Therefore, we can set v(r) = nyr!+¢/? for all r > 1 with n, = /2Ky + |a(0)|. 
Yang et al. [32] proved the existence and uniqueness of the strong solution 
of the scalar SDE (1) with Holder continuous diffusion coefficients presented 
in the following theorem. They proved the theorem based on [31, Yamada-— 
Watanabe theorem] and [10, Lemma 3.2]. 


Theorem 1. Let Assumptions 1 and 2 be satisfied. Then the SDE (1) with 
initial value x(0) = zp € R has a unique global solution x(t). Moreover, for 
all p > 0, there is a positive constant C that depends on T,, p, and xo such 
that 


a( sup |x(é)I?) <C, (12) 


O<é<T 


where E denotes the probability expectation under the probability measure 
P. 


For any real number R > |x(0)|, we define a stopping time 
Tr = inf{t > 0: |a(t)| > R}. (13) 


Lemma 1. Let Assumptions 1 and 2 hold. Fix any p € (0,00). Then, for 
any real number R > |(0)|, we have 


P(tr ST) < pp (14) 


where C’ stands for the generic positive real constant here in independent of 
R. 


Proof. By replacing € by tr A T in (12), we see 


la(tR AT)? <C. (15) 


Then, by Markov’s inequality, we have 


RPP(tR ST) < E(\2(tr) Plast) < E(le(tr AT)|) < C, (16) 


which completes the proof. 
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In the following sections, we consider numerical methods on a uniform 
mesh t, = nA forn=1,...,N with A=T/N for some N EN. 


3 An explicit two-stage truncated Runge-Kutta method 


In this section, we develop an explicit two-stage truncated Runge-Kutta 
scheme for the nonlinear SDE (1) with a superlinearly growing drift coef- 
ficient and a continuous Hélder diffusion coefficient. It is worth mentioning 
that in this regard, we adopt the idea of constructing the truncating functions 
of Mao [22] and Li, Mao, and Yin [18] to construct the new method. First, 
we outline some notations and preliminary results of the truncated methods, 
which will be used in the following sections. Further details can be found 
in the literature [21, 22, 11]. Let h : (0,1] + (0,00) be a strictly decreasing 
function such that for a constant h > 1 


A(1)>v(1), AM4a(A)<h, lim h(A) 


oo, for all A€ (0,1). (17) 


For example, we can consider h(A) = n,A~ with np, > nm, for any € € 
(0,1/4w), where w > 0. For a given step size A € (0,1], let ka: R-R 


z 
| ’ 


denote the truncation mapping defined by Ka(z) := (v~'(h(A)) A zl) for 
all z € R. We set 4 =O if z= 0. It can be easily deduced 


|z| 
lka(z1)| < |z1], |KA(Z1) — KA(Z2)| < 2/21 — z2|, ‘for all 21, 20 ER’. 


(18) 
Accordingly, we define the truncated coefficient by 


aa(z) := a(Ka(z)), (19) 
for all z € R. It is obvious from (3), (10) and (18) that 
laa (z)| < v(v7*(h(A))) = (A), for all z € R?, (20) 


and 
lan (21) — aa (z2)|? < 4K (1 + |zi/? + |22|")|21 — z2l?, (21) 


for all z1,z2 € R. 
Remark 2. From (10) and (18) it can be deduced for all z € Rand A € (0, 1] 
that 


laa (z)| < Jaa (z) — a(0)| + ]a(0)| < v(Y*(A(A)) Ka (2)] + |a(0)| 
< Ch(A)(1 + |zI). (22) 


Also, from (8) we can easily write 
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|b(z)|? < CC + |zl?). (23) 


In fact, relations (22) and (23) are the same inequalities (2.5) and (2.6) 
presented in Remark 2.1 in [33], which play a fundamental role in determining 
the optimal convergence rate of the partially truncated EM methods. 


The truncated function aa preserves the relation (9) in Remark 1 for all 
A € (0, 1]. We describe this in the following lemma. 


Lemma 2. Let Assumption 2 holds. Then for all p > 2, there is a constant 
K such that for all A € (0, 1] 


—l 
zaa(z) + lol)? < Kx(1+|z|?), for all zeR, (24) 


where Ka = $(Mi v [45] v [(e - 1)M3]). 
Proof. For any z € R with |z| < v~!(h(A)) from (7) and (8), we have 
zag (2) +P |b(2))? = 202) + P= *[o(2)P 
< (M1 V [(p — 1)M3]) (1 + |z/?). (25) 


Furthermore, for any z € R with |z| > v~+(h(A)), we can write 


Pee P= ole)? Say 


|2| 


=f (AG) a Lu 1(n(A)) a(vt(a(A)) = ) 


|2| |2| 


(1+ 1(R(A)))) + @— MBC + [ek*) 


a(y"(h(A)) 7) + (@— MRC + (2) 


)) + @= MBC + (2?) 


<3 (Mv [agay) Y le- Dg) a+ ee) 


(26) 


The inequalities (25) and (26) imply the required assertion (24) easily. 


Next, we construct our numerical algorithm TSRK2 to approximate the 
exact solution (1). For any given step size A € (0, 1], define 


Y (to) = xo 

Za(tn) = Ya (tn) + Ada, (Ya(tn)), 

Ya(tn+1) = Ya(tn) + A(araa(Ya(tn)) + a2aa(Za(tn))) + b(Za(tn)) AW, 
(27) 
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for n = 0,1,...,.N, where t, = nA and AW, := W(tn41) — W(tn). Here 
6 € Rand ay, a2 € [0,1] with a; +a2 = 1 are free parameters of the TSRK2 
procedure. As an example, for 0 = 0, we obtain the explicit truncated EM 
method [21]. Thus, the presented class of TSRK2 methods turns out to be 
a generalization of the truncated EM method. It is worth noting that when 
the truncated function a, is replaced by a, the TSRK2 method with a; = 0, 
az = 1, and @ = 0.5 is the SRK scheme DDIRDI1 [3] and reduces to the 
midpoint rule when applied to an ordinary differential equation. 

We now form a continuous-time version of the TSRK2 procedure (27). 
To do this, we first for any fixed step size A € (0, 1] set 


Ya(€) = D0 Ya(tn itn tn (6), Zal6) = do Zan ita tna), ‘for all € > 0. 
n=0 n=0 


é g 
uals) = ¥a(6)+ f (araa(Va(6)) +anaa(Za(c)))ac+ f 6(Za(o)aW (0), 


: : (28) 
or equivalently 


g € 
ale) = 20+ [(araa(Va()) + area(Za(c)))ac + f w(Za(Qaiv td): 
(29) 


3.1 Moment bound of the TSRK2 method 


In this subsection, we will state a new result showing the uniform bounded- 
ness of the solutions of the TSRK2 method (27). First, the following lemma 
shows how to conclude that the values Ya(€) and Za(€) are close to ya(&) 
with respect to the L?-norm. 


Lemma 3. For any p > 0 and any T > 0, there is a positive constant number 
Cy such that for every step size A € (0, 1] 


slya (€)—Ya(€)|PVElya (€)—ZalE)? < CpA?/?(A(A))?, for all € € on 
30 


Proof. We first prove the lemma for any p > 2. In this case, for a given 
€ € [0,T], there is an unique n > 0 such that tp, < € < tn41. By (28) and 
from the Hélder inequality, [20, Theorem 7.1], (8), and (20), one can conclude 


is 7 7 ig P 
Iva (6) ~ Yale)? <3°* (ad [Jaa (Va (6) Pa (31) 
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+ aA? 1 


g , oe F 
sf laa (Za(¢))|"do + A | \o(Za(C))|"a¢) 
<3?-1 (a A*(n(A))? + agA?(h(A)) 
+ 2?-! Ma AP/? (1 + |Z (tn) ?) 


Moreover, for every ¢ € [0,7], we can find a unique positive integer k such 
that th <€ < tg41. So, by (20) and (27), we have 


Za (C)-Ya()|? = 


Za (te)—Ya(te)I? = |O/PA? 


laa (Ya(te))|? < |aPAP(A(A))?. 

(32) 
On the other hand, by relations (8), (29), and (32) and the Doob martingale 
inequality, we have 


“Iya (€)|? <3?-* (aol? + 


g 2 g 
ff laa (Ya(¢))|Pd¢ + agE 
0 0 


Race sf b(Za(o))|Pac) < o(1+ (may)? [ 
7 € 
<o(1+(may?+ f 


2 € 
|Z (C) — Ya(C)|Pde + : A 


0 
(33) 
From (32) and (33), we can conclude 


yo 7 7 é 
sup El¥a(u)|? < sup Blya(u)|? < O(1+(h(A))?+ | 
O<u<é O<u<e 0 


By the Gronwall’s inequality, we can deduce 


sup E|¥a(u)|Pd¢). 
O<u<¢ 


sup E|¥a(g)|? < C(1+ (h(A))). 
0<€<T 


So, from (32), (34), and (17), we can write 


sup E|Za(é)|? < 2° *|6|PA?(h(A))?+2°-* sup 
O<esT 


‘|¥a(@)l? < C(1+(h(A))?). 
O<é<sT 


By inserting (35) in (31), we have 


"Iya (€)—Ya(E)? < C(AP(h(A))? + A7/? + AP/2(H(A))?) < CAP (A(A))?. 


(36) 
To prove the lemma for any p € (0,2), fix a number p > 2. Then by the 
Holder inequality, we can write 
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p/p 


sya (€) — Ya (6)? < (Blya(e) - Ya(eyl’)””” < (ca?2qa))?) 
= CpA?/*(h(A)). (37) 


So, from (32), (36), and (37) the proof of theorem is complete. 


Lemma 4. Let Assumption 1 and 2 hold. Then, we have 


sup sup Elya(é)|P<C, forall p> 0. (38) 
0<AK<10<E<T 


Proof. Let us fix A € (0,1] and € € [0,7]. Using the It6 formula from (29) 
for any p > 2, we can write 


é 
Iya (€)I? <lool? + f Plya (9)? *ya(¢)b(Za (6) dW (¢) 


g 
+ [ruler (ra(0 (a1aa(¥ (6) + e200 (Za(6))) 


+ Pt wzaior) ae, 
for all € € [0,7]. Since a, + a2 = 1, we have 
Iya (E)/? Slaol? + P1(€) + P2(g) + P3(€) 


g 
+ [ plya(Ol2ya (C)b(Za(Q)aw (0), (39) 


rie) = [rive (6P-* (ZalCaa(Zaley +> ozalonr)ac 
rag = [plo l-*(va(6) ~ Za(6)aa(Za(6)ae (40) 
rsa f * plua()P-2va(6)(aa(¥a()) - aa (Zale) ac 

By the facts that plya(¢)|?~*ya (¢)b(Za(Q)) is Fe measurable, we have 


S]ya(f)|? < Elao|? + Ei (€)) + E(P2(g)) + E(P3(€)). (41) 


Next, we try to estimate the values E(T,(€)), E((2(€)), and E(13(€)) in 
(41). By the relations (20), (24), and (27) and the special form of Young’s 
inequality (2), which reads 
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a 225 < ge + = 2) for all z1, z2 > 0, (42) 
Pp Pp 


we can approximate I'\(£) as follows: 


€ 
(11(g)) s sf ((p— 2)|y(O)P + 2KP/7(1 + [Za (Q)?)?/) de 


<(2K2)?/?7 + 230/2-1 KPT Ig PAP (A(A))? 


g € 
wo) | ‘Iya (C)[Pac + 287/2-Lice/? | s1¥a(Q)Pdc. (43) 


For T'2(€), with the help of Young’s inequality (42), relations (10) and (20) 
and Lemma 3, we obtain 


g T 
s(T2()) < (p— 2) | ‘lya(C)[Pa¢ + 2(A(A))?/2 / slya(¢) — Za(O)[P/2ac 


0 


é 
Si 2» | clya (C)IPaC + 27 Cy jo(h(A) PAP. (44) 


For ['3(€), the Young inequality as well as (10), (18), and (20) yields 


: p/2 
B(T3(@)) < an [((p—2)Iua(c)” + 2hva(C)/*Iaa(¥ (6) ~ aa(Za(e))|””) de 


€ - 
< a (p—1) | slya(C)IPd¢ + an (h(A))? [ Ina (¥(C)) — wa(Za(0))|PaC 


€ 
< a (p—1) | sJya(C)[Pdc + 2ay 4? T(h(A))2? A? 


§ R 
< (p- lay | sup Elya(u)|Pdu + 2|0|?a,Th?. (45) 
0 O<u<¢ 


Inserting (43)—(45) into (41), we have 


lua <c(a+ f * sup Elva (w)Pa¢). 


O<u<¢ 


Therefore, we can write 


€ 
sup we sc(i+ f sup Bly(u)|Pdc). 


O<useé O<us¢ 


It can be deduced from the Gronwall’s inequality that 


sup Ely(u)|? < C, 
0<u<é 
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for a positive constant number C' independent of A. For 0 < p < p, according 
to the Lyapunov inequality, (38) still holds, which completes the proof of the 
lemma. 


Remark 3. For every € € [0,7], there is a single integer n > 0 such that 
ty <&€ <tn4i1. From (17) and (27) we can deduce that 


‘|Z (€)|? =ElZa(tn)|? < 2°? (El¥a(tn))|? + A?|6|?E|aa (Ya (tn))|”) 


<2p-1 ( E|Ya (tn) |? + A3?/49 ne). 


Therefore, by Lemma 3 we can write 


sup sup E|Za(E)|? < C. (46) 
0<A<10<€<T 


In addition to (13), for any real number R > |x(0)|, we define two other 
stopping times 


og := inf{t >0:|Ya(t)|>R}, 0 :=inf{e>0:|Za(t)| > RB}. (47) 


Lemma 5. Let Assumption 1 and 2 hold. Fix any p € (0,co). Then, for 
any real number R > |a(0)|, we have 


PO <T)VPOY <T)< = 


a) d 


where C stands for the generic positive real constant here in independent of 
R. 


Proof. The proof of this lemma is similar to that of Lemma 1. Namely, 
replacing € by 04) AT in (38), we see 


ya(Oe AT) <C. 


Then, by Markov’s inequality we have 


RPP(OR? S$ T) < E(lya(OR’)PLycp) S Eya(g’ ATP) < 


Moreover, it follows from Remark 3 that E |Z (08 dA T)|? < C. Therefore, 
as above 


RPP(OR? ST) SE(IZaOe I Lycp) S E(Za (Oe AT?) < 


which completes the proof. 
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3.2 Convergence of the new method at time T 


In this section, we study the effectiveness of the method TASK2 in solving 
problem (1) and obtain the corresponding convergence results at finite time 
T. Using the method of Yamada and Watanabe, for each 6 € [1,co) and 
€ € (0,00), we choose a nonnegative continuous function ys. : [0,co) > R 
such that [15, 31]: 


© x5e(E) < TOE if € € [e/d,€],  x5e(€) =O otherwise; 


be Sos NoelE)GE =1. 


Moreover, we approximate the function € :—> |&| by the function ws, defined 
by 


\El pCa 
mge(€) = | i: VXs5e(C2)dGad¢,, forall € ER. (48) 
0 Jo 


Below we outline some properties of the ws, function that will be used in the 
following sections. Further details can be found in the literature [32, 33]. 


Lemma 6. Let 6 € [1,00) and let € € (0,00). Then for all€ ER 
1. IE S w5e(€) + €, 


3. §-(€) = Xse(IE|) < Wins le/6<lel<e}) 
4, Zi > 0, for all € ER \ {0}. 


Remark 4. Since supgeR|™5.(§)| < 1, under Assumption 2, we can easily 
conclude 
W§e(21 — 22)(@(z1) — a(z2)) S$ Ailes — al, (49) 


for all z1, z2 € R; see [32]. 


Define the stopping time 6a rp := TR A gt) A g®), where Tr and g® 
for i = 1,2, are defined by (13) and (47), respectively. Then, the moment 
deviation between x(€ A Ba.r) and ya(& A Ga,p) is estimated as follows. 


Lemma 7. Consider the initial problem (1), which satisfies Assumptions 
1 and 2. Suppose that R > |xo| is a real number and that A € (0,1] is 
sufficiently small such that v‘~!)(h(A)) > R. Then there exists a constant 
C such that 


| 
5 
B BR 
+ 
b 
él 
= 
cs 
We 


C if a=0 
ealE Pal < | ), fa=0, (50) 


C(Azh(A))™, if 0<a<0.5, 


where 
ea(é) :=2(€)—-yalé),  for0<€<T. 
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Proof. We note for0 <¢ <€ABaR gars [¥at )| 
Therefore, we conclude by the condition v(-)) (h(A)) 
of the truncation function a, in (19) that 


aa(Ya(¢)) =a(¥a(¢)) and — aa(Za(G)) = a(Za(¢)). 


It therefore follows from the It6 formula that 


R and |Za(¢)| < BR. 


< 
> R and the definition 


flea (EA Ba,r)| Se + E(wse(ea(E A Ba,r))) (51) 
EABA,R 
ae +B [8 whclea()[ale(@)) ~ ara(¥a(@)) ~ analZa(O)] ae 


1 
2 


+ 


EABA,R 
: i, of. (ea(C)) [ble(¢)) — (Za (0))]2a¢ -= + Th + M12. 


In what follows, we attempt to estimate the values II, and Hg, in (51). Since 
Qa, + a2 = 1, by re-arranging we get that 


EABA,R 
Il, <E i: mhe(ea(6))[a((C)) — aya (¢))] aC 


EABA,R 
+E | lari. (ea(C))||a(ya(C)) — a(¥a(C)) dc 


EABA,R 
+ aE i: lo4-(ea(C))||a(¥a(0)) — a(Za(Q) ac 


By Assumption 1, Lemma 6, and Remark 4 we can write 


EABA,R 
Tl, <H,E i lea(Q))|ac 


EABa,R 
tal ee ff (1+ lya(Q)l? + 1¥a(¢ IP)? |ya(¢)) — Ya(¢)|d¢ 
EABa,R 
saa Ie ff (1+[¥a(Ol? +1Za(Ql)2|¥a(O) — Za(Qlac. 
0 


So, by the Hélder inequality, Lemmas 3 and 4, we have 


EABA,R ‘ 
Th < HE | lea(C))|d¢ + CA4A(A). (52) 


As for Ig, it follows from Assumption 2 and Lemmas 3 Lemma 6 that 


EABA,R 1428 1 F 
< HSE —Z —_____ | 
Ty < Ay | (a(¢) — Za(6)) lea(@Ind{e/8Slea WISe ¢ 


PoHeeog gapyas pé 1420 
! -Z 
In(6) eln(6) | Iya a(¢)| 


d¢ 
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1+2a 


< 2° 13T| “ A?h(A)) (53) 


Substituting (52) and (53) in (51) and by the Gronwall inequality, we have 


2a 


€ 6 1 1420 
in(ay + cinta (47 (A) ). 


Hlea(E A Ba,n)| <C(e+ A?K(A) 4 


In the case when a = 0, we set 6 = A~# ande= — Tray which implies 


lea(E A Ba,n)| $ C(— pray + A4M(A))4). (54) 


If we take 6 = 2 and e = Ath(A) for the case where a € (0,0.5), then we 
have 


slea(E A Ba,r)| < O(A#A(A))™. (55) 
The inequalities (54) and (55) prove the desired. 


Theorem 2. Let conditions in Assumptions 1 and 2 be fulfilled and let p > 1. 
Let Ra := (Azn(A)) P79 for any A € (0,1). If there is a positive real 
number A* € (0, 1] such that 


v-l(h(A)) > Ra, for all A€ (0, A*]. (56) 


Then, there exists a positive constant C' independent of A such that 


C(-piy + At(h(A))?), ifa=0, 
lea (T)| < ( a - ( ( )) ) nwa (57) 
C(Azh(A))™, if a € (0, 5). 
Proof. We first divide the left side of (57) as below: 
lea (T)| = Elea(T) 146, 2>7}| + Elea(T)1t6.,2<7}I (58) 
By Young’s inequality (2), Theorem 1, and Lemmas 4 and 5, we obtain 
A2h(A), pa 
lea(T) Lie, n<T}I < en) tlea (T)|? aC 1 1/(p—1) P(Ba R S T) 
e p(Azh(A)) 
< CA?h(A). (59) 


Since v~1(h(A)) > Ra, substituting (59) in (58) by Lemma 7, we obtain the 
result of the theorem. 


Remark 5. This theorem shows that the TSRK2 method has an order of 
convergence close to a for a € (0,0.5). This is theoretically almost opti- 
mal if we recall that the classical EM method has a convergence order of 
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a. However, the condition (56) could sometimes make the TSRK2 method 
impractical; see [33] for more details. 


In what follows, we use the mathematical techniques developed by Yang 
and Huang [33] to remove the imposed condition (56). In this context, 
we replace the condition A'/4h(A) < h by the more general condition 
Al/2h(A) < h for all A € [0,1]. Without limiting generality, for simplic- 
ity, we use here h(A) = n,A~? with np, > 0. 

Lemma 8. Let Assumptions 1 and 2 be satisfied. If h(A) = n,A~? with 
Nn > 0, then for every A € (0,1] and every p > 0, it holds that 


sup E| sup |ya(é)l"] $C, (60) 
0<A<1 ‘o<é<T 
and 
( sup Elya(é)— Ya()I’) V (sup _Elya(€)— Za(@)?) <CA®, (61) 
O<eé<T O<€<T 


for all real positive number T. 


Proof. We fix A € (0, 1] and use the same notation as in the proof of Lemma 
4. By using the It6 formula, we can write that for any p > 2, 


sup, valOI? <laol? +P (L) + Pal) +1s(7) 
é 

+ sup | [riua(OPva(c)oZa(c)aw(o)}. (62) 
0<é<T'/0 


So, by the Burkholder—Davis-Gundy inequality and the linear growth condi- 
tion property of b in (8), we have 


r( sup |ya(€)l) <leol” + Ta(T) + P2(7) +P (7) (63) 
O<é<T 


ange o([ sup _lva lr? [ware + \Za(€))?) df] ia 


In what follows, we estimate the values E(T'\(7)), E({2(T)), and E(T3(T)) 
in (63) with the more general condition h(A) = 7,A~2. For E(1,(€)), we 
obtain from (43) 


T Te. 
rig se(i+ f° BlyaoPac+ f° Elvacyrac). 4) 


As for T'2(€), we obtain by Young’s inequality (42) that 
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T T . 
(Pa(6)) <2) f Blya(c)Pa +28 | Iya(G) ~ Za(6)] laa ¥a (6) Ed 
0 0 
T 
<CA-* f (Elya(6) — Za(6)|P)? (1 + El¥a(C)I?)26 


0 


T 
+(p—2) i ‘Iya (OPAC. (65) 


Since (h(A))??A? < 7? h?, following an approach very similar to that used 
in (45), we can show that 


T 
“(T3(€) < C(1 + | Elya(C)|Pdu). (66) 


For any ¢ € [0,7], there is a unique integer n > 0 such that tn < ¢ < tn41. 
Then we can write 


|Za(6) — Ya(¢)|? = E|Za (tn) — Ya(tn)|? = A?[O|"E|aa (Ya (tn))|? 


and 
Slya(¢) — Yad)? =Elya(¢) — Ya(tn)|? < CAP (mAs [laa (Ya (tn))?] 


+ a2AFE[laa(Za(tn))?] +E[I(Za(tn))!"])- (68) 


Therefore, by Remark 2 and relations (67) and (68) we can obtain 


ya(6) — Ya(¢)? < CA? (1+ E|Ya(tn)|? + E|Za(tn)/”) 
< CA? (1+E|¥a(¢)|?). (69) 


x 
vy 


From (67) and (69) we can simply conclude 


lya (6) — Za(C)/P < CA2 (1+ E|¥a(¢)I?). (70) 


Substituting (64)—-(66) and (70) into (63) and applying Young’s inequality 
that is 


Zyzq <= + —, for all 21,22 € R, forall 6>0, (71) 


with 6 = 8pMb, we can write 


r( sup |ya(é)l?) <C(1+ | * Blya(G)lPac + i * EIYa(6)IPaC) 


O<é<T 
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+ 5E( sup |ya(€)’) 


O<é<T 


+ 32p° MZE( [src + |Za(€)/?)ag) 


< r( sup |ya 6") 


ro(it f Biaiorac+ f° Biva(cra), 


which implies 


( sup Iya (€)I”) = ot+f o[ sup Iva(s)"Jds) 


0<e<T 0<s< 
Then, the Gronwall inequality implies the relation (60) for p > 2. For p € 
(0, 2], (60) still holds as desired due to the Lyapunov inequality. Moreover, 
(61) is obtained directly by substituting (60) into (69) and (70) , which 
completes the proof. 


Theorem 3. Consider the initial problem (1) that satisfies Assumptions 1 
and 2. Let v(r) = mr!*? and h(A) = n,A~? in which p > £. Then there 
exists a constant C such that 


C(-aty At), if a =0, 
lea(T)|<4~\ mA TAT” O (72) 
CA®, if 0<a<05. 
Proof. We split the left side of (72) into two parts 
lea (T)| = Elea(T) Iya, n<7}| + EleaT)6.,2>7}1- (73) 


Using the Young inequality, we have 


A2 aug pad 
lea (T)I SS Blethen 
lea( Migansti] S 545 jea(T)| 54S 


A- 7070 P(Ba.p <T). (74) 


For any R > |xo|, by Lemmas 1 and 5, we can write 
RPP (Ba < T) <C, 


which follows that C 
P(Bar<T)< Ree (75) 


If we substitute (75) into (74) and use Theorem 1 and Lemma 8, we get 


bd 
A 2@FD ) 


7 i 
Hlea(T)Lre,. p<T}I < c(a} + ~ Retz 
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By choosing R = v~1(h(A)) = (2A72) 1, we have 


blea(T) Ita, n<T}| < CA?. (76) 


In what follows, we try to estimate the second terms in (73). Since R = 
v—l(h(A)), if 0 < ¢ < Bar, then |¥Ya(d)| V |Za(Q)| < v71(h(A)) for all 
A € [0,1]. By the definition of the truncated function (19), aa({Ya(¢)|) = 
a(|Ya(¢)|) and aa(|Za(¢)|) = a(|Za(¢)|). Using the same notation as in the 
proof of Lemma 7, according to the It6 formula we get 


flea (E A Ba.r)| <e4+I, +12, (77) 


where II, and [2 are defined in (51). According to Assumption 1, Lemma 6, 
and Remark 4 as well as the Holder inequality the following equation holds 
for TI, 


€ 
Th, <Hy i ‘lea (CA Ba.n))|aC 
(0) 


ya(6)) — Ya(Q)|?)2a¢ 


EABA,R 1 
+ vm [ (1+ Elva? + El¥a(Ol*)? E 


EABA.R 1 1 
+o2/Ki. sf (1+ E|¥a(0)l? + ElZa (O°)? (E|¥a (6) — Za(6)|")2 de. 


From Theorem | and Lemma 8 we can thus conclude 


§ a 
Tl, < of slea(CA Ba.p))\de + C3. (78) 
0 


As for IIg, it follows from (6) and Lemmas 6 and 8 that 


' EABA,R oe 1 P 
< D — —-— 
Mla < HBE [  (0(6) - Za()* ep ltcecteaconcat 


haa al cae Oe 1420 
may t tary E ff lua) - Zatol ag 
E20 bA1/2+e 


< CD - ney)” (79) 


Substituting (78) and (79) into (77) and applying Gronwall’s inequality yield 


e+ e* §AV2to 1/2 
ne) Ray 


lea (EA Ba,r)| < C( 


In what follows we divide the proof into two cases a = 0 and a € (0, s). In 


the case when a = 0, we set d = A~4 and ¢ = 1, implying 


slea(E A Ba,p)| < C-ray ses, (80) 
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In the case where a € (0,5), we set A = 2 and e = A2, which implies 


slea (EA Ba,r)| < CA®. (81) 


From (77), (80), and (81) the proof follows directly. 


4 Convergence rate over a finite interval 


Sometimes we need to approximate path-dependent quantities, for example, 
the value of the European barrier option value. In these cases, we need a 
stronger convergence result such as 


Jim E( sup |e(@) ~ valé)) = 0 


Lemma 9. Consider the initial problem (1), which satisfies Assumptions 
1 and 2. Suppose that R > |x| is a real number and that A € (0, 1] is 
sufficiently small such that v(-)(h(A)) > R. Then for all 0 < € < Tia 
constant C’ exists such that 


e( sup |a(6) — ya(6)I 


0<C¢<E ~ : Ace 


C(AtA(A))** , if a € (0,4). 


Proof. By the It6 formula, we can write 


lea(E A Ba,r)| Se + (@5,<(ea (EA Ba,r))) 


where, for € > 0 


EABa,R 
eKO= wh (ea(0)(b(«(0)) —B(Za(O))aW(d). (83) 


Since v(~)(h(A)) > R, if 0 < ¢ < €A Bar, we have |Ya()| V |Za(0)| < 
v\-Y(h(A)), which gives aa(Ya(¢)) = a(¥a(¢)) and aa(Za(¢)) = a(Za(¢)). 
If we use the same technique as in the proof of Lemma 7 and apply Assump- 
tion 1, then we have 


EABA,R 
lea(€ A Ba,a)| <e+ Hh | lea (Q)lae (84) 
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EABA,R 
+Vm | (1+ ya? + Ya OQ!)?|yaQ) — YaOlac 


EABA,R 
+avKi | (1+ |Ya(QP + Za)? |¥a(O) — Zale 
Vad 6 Ota 
“  In(6) 
2a 772 EABA,R 
ama | yale) =a) ae oes ae: 


" eln(6) 


By the Holder’s inequality and Lemmas 3 and 4, we have 


:( any p leah Ba, n)l) <e+ Ay i ( sup lea(wA Ba,n)|)dG + CA2H(A) 


O<u O<u<¢ 
920 F220 920 F725 : € 1420 
itd) Seng) Ef bal Zalo*a 
=a *( sup |$5..a()I): (85) 


Next, we try to estimate the value E (suppcucelSs<.a(w)) in (85). By the 
Burkholder—Davis—Gundy inequality and Lemma 6, we can write 


(sup ISs-a(u)l) <4v3E( f ee ohe(ea (6))F|(2(G)) — 8(Za(C)) Pac) 


O<usé 


<xettinae( ["o()~ 2alonae)’ 


5 EABA,R A 
+2ettHne([°  ea(6)- Za(c)[*2d0) 
0 
In fact, by Assumption 2 we have 


|b(a(¢)) — b(Za(¢))|” < Hale(¢) — Za(¢)|2 
< 2° Ha|x(¢) — xa(¢)|'+?* + 2° Aa|axa(C) — Za(Q)|*7*. 


Therefore, by Lemma 3 we have 


:( sup |55,,A(u) ) <c (At Va(d)) 


O<u<é 


1+2a 


+otine( f NN" (6) — 2a (O22aC)* (80) 


In the rest, we divide the proof into two cases. 
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Case 1: Suppose that a = 0. In this case, according to Lemma 7 and from 
the relation (86), it follows that 


7 ; + 1 z 4 3 
( sup ,!S.(%)]) <o(a Vi) + (RH aay + PHA) ) ) (87) 


Substituting (87) into (85) and choosing 6 = A~# and ¢ = —1/In(A), we get 


g 
x( sup lea(wABa.a))) <Ah [°B( sup lea(u A Ba,n)l)ac 
O<usé 0 O<us¢ 


-c(ainay + A2(h(A))? : 


which implies 


g 
B( sup lea(wABa,a)l) <A [B( sup lea(u A Ba,n)l)ac 
O<useé 0 \0<u<¢ 


1 1 1 3 
——— + A7(h(A))? ]. 
+ (aay + Atay?!) 

Case 2: Suppose a € (0, s). In this case from (86) we have 


( sup [55,<,a(u)l) <c(a 
O<usé 


(88) 


Es (A)) 142a 


(28 leal€A 8,2) 2 ( [ ee bea (Pd) ‘) 


By Young’s inequality and Lemma 7 we have 


(sup |S5c.a(w)l) <o((a* Vay) = ivanay*) 
O<uxé 


5B ( sup lealé ABa.n)l)). (89) 


Inserting (89) in (85) and by choosing ¢ = A2h(A) and 6 = 2, we get 


g 2 
:( sup lea (uABa,n)l) S 2m, | ( sup lea (wABa,n)|) a¢+C (VAa(A)) : 
O<u<é ) Osus¢ 


(90) 
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Hence, by using the Gronwall inequality on (88) and (90) we have the required 
assertion. 


Theorem 4. Let v(r) = mr't? and h(A) = m,A~“* in which p > §, and 


gull 3 
€€ (0,min{ 7. on 
for any w € (0,2). Let Assumptions 1 and 2 hold, and also let A* be suffi- 


ciently small such that Aa < —4_. for all A € (0, A*]. Then 


In At 

1 2 : = 
-( sup J2(6)-vatol) <4 Gms) ae. oi 
0<¢<T C Ara*(1-20) 4 ATs), if #& (0,2), 


Proof. For any R > |x|, we can decompose the left side of (91) into two 
parts as below: 


E( sup |x(6)—ya(6)l) <E( sup Jea(wABa,n)l)+E( sup lea (6)IItaa,x<7}): 
O<¢<T O<usé O<¢<T 
(92) 


Let p> 1+ 322 be fixed. By the Young inequality, that is 


2we 


Op, P-1 _ if(p-1) 


2122 < ov ++ oat eat) 2 , for all 21,22 € [0,+00), for alld > 0, 
Pp 


and Lemmas | and 5, we can deduce 


) p-1 
: =K Pp eee < 
(sup Jea(Olltes.xsr) SS (sup eal) ) a poipai (Ba.r ey 
—1 
gO) a. 2 (93) 


~ p + 61/P-1 Rp’ 


By choosing 
6= AG=ar. R= A=U=9oPe—)) 


we see that v(—(h(A)) > R. Therefore, by substituting (93) into (92) and 
applying Lemma 9, we have 


1—2we 


C((adsy tA )F +A"), if a=0, 


C( A207 A206) +A‘), if a€(0,4 


E( sup |2(¢)—ya(6)| 


0<¢<T 


sr 


| 


which completes the proof of the theorem. 
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log,, plot 


I-3—TSRK2 method 
|-- Truncated EM method 
|} Partially truncated EM method 


0.67} ___ Reference curve: 0.5log,, (-1/In()) 


E(supy <ecql *(€) - Yq (6)I) 


Figure 1: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 1 with « = 0.0001, w = 0.5, and a = 0. 


5 Numerical experiments 


In Section 3, we derived TSRK2 methods with free parameters a1, @2, and 
6. In this section, we follow [17] and set them to below: 

1 Wy 

ay =a, == and g=—, (94) 

2 Wo 
in which wo = 1.9 and w; = 0.8184. We confirm the performance of the new 
method in terms of accuracy and stability in contrast to the truncated EM 
method [21, 32] and the modified partially truncated EM method [33]. In 
this context, we investigate 


2000 


» (95) 


1 - . 
; = — (t,) — yt 
(sup let) — yall) B00 Da BRAY yl En) — VS) 


for a given A = T/N to measure the accuracy of the methods. 


Example 1. Consider the following SDE with Hélder continuous diffusion 
coefficient 


dex(é) = (Ar2(é) — Aaa¥()) dé + wle(Q)l**°aw(e), €€10,T], (96) 
x(0) = 2 ER, 
where \, € R and Ao, p € R*. 


The drift and diffusion coefficients are a(z) = Az — Azz? and b(z) = 
plz|2te, respectively. We can easily show that, for 21, z2 € R, 
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log,, plot 


[—=—TSRK2 method 
|-©--Truncated EM method 
~~ Partially truncated EM method 
|----Reference line, slope=0.5 


‘oe 


0.5 - 


-0.5,- 


E(supy ecql *(€) - Yq (6)1) 


Figure 2: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 1 with « = 0.0001, w = 0.5, and a = 0.25. 


|a(z1) — a(z2)|? < Ki(1 + laal* + |z2|*)|aa — 221”, 


where Ky := 3max{A7, 243}. In other words, Assumption 1 is satisfied with 
p =4. Moreover, we have 


(z1 — z2)(a(z1) — a(z2)) < Ar(z1 — 22)", [b(21) — b(z2)| < pla — 22? *°. 
That is to say, Assumption 2 is fulfilled. Concerning (10), we set n, = 4K}. 
Then, we have 


2 


Al Zz — AalzZ 
re lata) — atza)l <5 ( oy < mr 


0<|z1|V|za|<r |21 = Z| 


and 

sup |a(z)| < (|Ar| + A2)r? < mr®, 

lz|Sr 
for all r > 1. To apply Theorem 4, we set v(r) = n,r? and h(A) = n,A~“.. 
We can therefore conclude that the truncated Runge-Kutta solution (27) 
satisfies 


1 


C 1 2 if a=0 

( s - < {Cl m=)" 

(sup Je —walQl) $4 Ra ay 
be cl a2 +A3), if ae(0,5). 


Let us now test the efficiency of the TSRK2 method for problem (96) in 
comparison with the truncated EM method [21] and the modified partially 
truncated EM method [33] for different values of the Hélder parameter a. 
Since we do not know the exact solution of (96), we search for a numerical 
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log,, plot 


—=—TSRK2 method 
--©-Truncated EM method 
me Partially truncated EM method 
----Reference line, slope=0.5 
0.5- 


“oe 


-0.5,- 


E(supy <ecql X(€) - Yq (6)I) 


Figure 3: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 1 with « = 0.0001, w = 0.5, and a = 0.45. 


solution with A = 2~!§ using the truncated EM method and use it instead 
of the exact solution. Figures 1-3 show simulations of (95) for the TSRK2 
method with parameters (94), the truncated EM method, and the modified 
partially truncated EM method with w = 0.5. Here we take Ay = —1.5, 
A2 = 10, and w = 1.5 as parameters of (96) and apply the methods over 
0 < € < 2 with step sizes, A = 2'~? for i = 1,2,...,9. In these figures, the 
error of the methods is plotted as a function of seven-step sizes A in logio 
for different values of a (0,0.25 and 0.45, respectively). 

Figure 1 shows that, compared to the reference curve 0.5 logio(—jazay) 
the TSRK2 method gives better convergence results than the truncated EM 
method and the modified partially truncated EM method for step size 2—° < 
A < 2-!. However, for sufficiently small A (e.g., A € (0,27°}), all methods 
have the same slope values as in Theorem 4. Increasing the value of the 
Holder parameter a, it can be seen in Figures 2 and 3 that for the step size 
2-5 < A < 27!, the TSRK2 method has a better convergence rate than 
the other methods compared to the reference line with slope 1/2, but for 
sufficiently small A < 2~®, all methods have the same convergence rate close 
to 0.5 as in Theorem 4. The simulation results clearly show that the TSRK2 
method confirms the convergence results stated in Theorem 4 and that the 
new method is the most efficient method in terms of accuracy compared to 
the truncated EM method and the modified partially truncated EM method. 


Example 2. Consider the scalar nonlinear Its SDE with a one-dimensional 
Wiener process 

de(t) = (Ara(E) — A20°(€)) dé + nle(@)|2 aw), t € [0,7], 
x(0) =1 


) 
)=1, (97) 
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—S—TSRK2 method Pe 
--©-Truncated EM method 
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1.5 


ab 


0.5 - 


E(supp <ecql *(€) - Yq (6)I) 


xe 
@ 


Figure 4: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 2 with « = 0.1, w= 0.5, and a=0. 


where A; € Rand Ag, € Rt. 


The drift and diffusion coefficients are a(z) = \1z — A2z° and b(z) = 
plz|2 te, respectively. We can easily show that, for 21, z2 € R, 


|a(z1) — a(z2)|? < Ko(1 + Jaa|® + |z2l®)[en — zal, 


where Ky := 3max{)?, 223}. In other words, Assumption 1 is satisfied with 


p = 8. Moreover, we have 
(21 — 22)(a(z1) — a(za)) < Aa(z1 — 22), |b(z1) — B(zo)| < plea — 20|2 4°. 


Concerning (10), we set y, = 4V/ Ko. Then, we conclude 


|a(z1) — a(z2)| Z ale +r 


5 ys mr’, 


sup 
0<|z1|V|zal<r |z1 — 22| 


and 

oe Ja(z)| < (lAa| + A2)r? < mar’, 

z\<r 
for all r > 1. To apply Theorem 4, we set v(r) = nyr° and h(A) =n, A“. 
Figures 4-6 show the simulation of (95) of the TSRK2 method for equa- 
tion (97) as a function of step size A for Holder a = 0, a = 0.25, and 
a = 0.45, respectively. Since there is no explicit solution, we use the trun- 
cated EM solution with A = 2—'8 as a suitable approximation to the exact 
solution. Here we take \; = —15, Ag = 20, and w = 3 as stiff parame- 
ters of (97) and apply the methods over 0 < € < 2 with different step sizes 
Cay eee ee 
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Figure 5: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 2 with « = 0.1, w = 0.5 and a = 0.25. 


Figure 4, corresponding to a = 0, shows that the TSRK2 method gives 
better convergence results than the truncated EM method and the modified 
partially truncated EM method, especially for the step size 2~° < A < 273, 
compared to the reference curve 0.5 log; (— may)! For the step size A = 2-3, 
for example, the new method has an error of 0.8655, while the truncated EM 
method and the modified partially truncated EM method reach an error of 
4.9260 and 25.1876, respectively. It should be emphasized that the modified 
partially truncated method EM is not applicable for A > 2~?. However, for 
sufficiently small A (e.g., A < 2~°), all methods have the same slope values 
as in Theorem 4. 

In Figures 4 and 5, we present the supremum norm of the methods for 
larger values of a = 0.25 and a = 0.45, respectively. These figures also show 
that the new method has a better convergence rate compared to the other 
methods mentioned, especially when the step sizes are not very small. 


6 Conclusion 


In this work, we have developed an explicit TSRK2 scheme for nonlinear 
one-dimensional It6-SDEs with one-sided Lipschitz and local Lipschitz drift 
conditions and continuous Hélder diffusion condition. We proved the mo- 
ment boundedness and convergence properties of the approximate solutions 
at a time T and in a finite time interval [0,7], respectively. We have demon- 
strated the efficiency of the new method using some numerical examples with 
different Hdlder parameters a € [0, 5) and compared our method with the 
truncated EM method and the modified partially truncated EM method in 
terms of their strong convergence performance. The numerical simulations 
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Figure 6: The supremum convergence norm (95) in [0,2] as a function of step size A for 
Example 2 with « = 0.1, w = 0.5 and a = 0.45. 


we performed confirmed the theoretical convergence results and also showed 
that the TSRK2 method is robust to changes in the Hélder parameters of 
the SDEs. 
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